Localized Entanglement in one-dimensional Anderson model 
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The entanglement in one-dimensional Anderson model is studied. The pairwise entanglement has 
a direct relation to the localization length and is reduced by disorder. Entanglement distribution dis- 
plays the entanglement localization. The pairwise entanglements around localization center exhibit 
a maximum as the disorder strength increases. The dynamics of entanglement is also investigated. 
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Entanglement is a kind of nonlocal correlation that 
only exists in quantum systems. Recent studies on en- 
tanglement are motivated by its potential applications in 
quantum computation, quantum teleportation and quan- 
tum communication. As spin system is a perfect one 
to realize quantum computer, many efforts focus on the 
entanglement in the Heisenberg spin model 0,0, Ising 
model in a transverse magnetic field 0, and itiner- 
ant fermionic systems yj. In quantum computer, one 
needs to control and measure on individual qubits. How- 
ever, in many possible physical implementations of quan- 
tum computer, the interaction between spin and spin (or 
qubit-qubit) is inevitable, and excitation or paticles can 
hop from one site to other sites. So it is hard to operate 
on one single qubit. To overcome this difficulty one ex- 
citation should be pinned on a certain site. In fact, it is 
well known that the localization can pin the excitation. 

On the other hand, being a fundamental concept of 
quantum theory, entanglement must be involved in many 
fields of physics. It is shown that entanglement is a indi- 
cator of quantum phase transition 0, 013 • The relation 
between entanglement and chaos, and the relation be- 
tween entanglement and localization are also discussed, 
and it was found that strong localization decreases en- 
tanglement [jJ.The quantum entanglement in condensed 
matter system, bose system and fermi system, and their 
connection with long-range order and spontaneous sym- 
metry breaking are also discussed |Ullll| . 

The effect of incommensurate coupling strength on 
pairwise entanglement was studied by considering the 
one-particle states of the Harper model ^| an d Frenkel- 
Kontorova model , respectively. The behaviors of en- 
tanglement during the phase transition and relation be- 
tween pairwise entanglement and state localization were 
revealed. Moreover, there is also a simple quantitative re- 
lation between the bipartite entanglement and the state 
localization 1141. 



Disorder is a common factor which exists in a large 
area of physical world. It is well-known that disorder can 
lead to localization from the Anderson's early work [lfij . 
which influences many properties of physical system such 
as the electric insulator. In this paper, we study the 
effects of disorder on entanglement of one-particle states 
in the one-dimensional Anderson model. We show that 
the localization decreases the entanglement sharing in 
one particle state and entanglement is shown to be an 
indicator of the localization in one-dimensional systems. 

In general, one-particle state belongs to a subspace of 
the 2 Ar -dimensional Hilbert space of N qubits. This sub- 
space is TV-dimensional and spanned by states with only 
one excitation. One particle state can be written as 
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this state can be realized in many quantum systems such 
as those of one spinless fermion or boson hopping on a 
substrate, and one magnon in Heisenberg spin systems. 

For a pure state of bipartite system, composing of 
two subsystems A and B, the bipartite entanglement 
can be measured by Von Neumann entropy, linear en- 
tropy, or other entropies. For mixed state of two qubits 
p, Wootters et.al |ljj had found that the entanglement 
of formation is monotonic function of its concurrence 
which is defined by C(p) — max(0, Ai — A2 — A3 — A4) 
where Ai,A2,A3,A4 are the square roots of the matrix 
P= p(o-y (g) a y )p*(o- y ® ay). For the general one-particle 
state, the entanglement between a pair of qubits, qubit i 
and j, can be quantified by the concurrence given by [l?! ] 
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One electron hopping on substrate potential can be de- 
scribed by a general Hamilton H = ^ + V. If we are not 
interested in the behavior of the wave function on length 
scales smaller than a lattice constant, the model can be 
described by the discrete Schrodinger equationwith with 
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a fixed lattice constant a 

h 



2ma' 



(ipi+i - 2tpi + ipi-i) + Vnpi = Eitpi, (3) 



which can be written in more comprehensive way as 

-t(ip i+ i + i/'j-i) + s l ip i = Eifc. (4) 

In the second quantized picture, the Hamiltonian can be 
written further as the following 
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where t is the nearest-neighbor hopping integral, mea- 
suring the probability for electron transfer from n-th 
site to its nearest-neighbor sites, it is chosen to be 
unit throughout this paper. Vi is the on-site poten- 



tial. 



and Ci are the creation and annihilation oper- 



ators of i-th local fermionic modes, satisfying the canon- 



ical anti-commutation relation 



Cj]+ = 0, 

,j . The general state of the electron hopping in the 
one-dimensional lattice can be viewed as a multiqubit 
state in the occupation-number basis, and thus the 
entanglement between two local fermionic modes can be 
discussed [l^l-In fact, the single electron model is equiv- 
alent to one magnon state of the XY model which is 
described by the Hamiltonian H = —tYln=i( a n a n+i + 
a n a n+i) + ^2n=i hnCn- So the entanglement properties 
of these two models are the same, we can also under- 
stand the effect of localization on the entanglement of 
pair spins (qubits). 

Consider the motion of electron in the one-dimensional 
Anderson model, the on-site potential can written as 
Vi = Vq + Xei, where e, is a Gaussian random variable, 



satisfying 



■>13, 



(e») = 0, and A is the disorder 



strength. To study global entanglement of the system 
and reveal the relation between entanglement and state 
localization, we use the entanglement measure given by 
the average concurrence |l2j 
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with M = N(N -l)/2. 

As mentioned in the introduction, the localization 
caused by disorder is typical and important in condensed 
matter physics. The direct result of disorder is the local- 
ization of state, i.e., the wave function of system, which 
exhibits a exponential spatial decrease |l9l | , 



ipi ~ ip io exp 



(7) 



where £ is the localization length, Xi is the site coordinate 
of wave function, and Xi is the center site of localization. 



FIG. 1: Average concurrence as a function of disorder 
strength parameter A. Solid line is the fit curve of exponen- 
tial decay. The system size JV = 1600 and the concurrence is 
averaged over 50 disorder ensembles. 



In the continuous limit, the average concurrence © 
can be written as an integral. Considering the fast expo- 
nential decrease of wave function of localized state, the 
average concurrence can be estimated as 
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where A is the constant. It is evident that the average 
pairwisc entanglement has a direct relation to the local- 
ization length in the one-dimensional disorder system. 
The localization length is the most important indicator 
of localization especially for the disorder system. When 
disorder becomes stronger, the localized length become 
small, the state of system is more localized and less en- 
tangled due to the above analytic result JSJ. 

To show the above result in more detail, we study 
the entanglement of ground state of the present one- 
dimensional disorder system by numerical calculations. 
In Fig. 1, the behavior of average concurrence of the 
ground state of Hamiltonian(5) versus disorder strength 
is presented. When the disorder is absent, i.e., A = 0.0, 
the average concurrence exhibits a maximal value. The 
increase of disorder strength leads to the decrease of av- 
erage concurrence. It is shown that the decrease of con- 
currence fits the second exponential order 
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with Bx = 0.00872, D 1 = 0.10611,52 = 0.00139, and 
D-2 = 1.2166. This numerical result indicates that 
the disorder has great effects on quantum entanglement, 
agreeing with the analytic result JSJ.In the extended case 
of Hamiltonian(5) , the system has a maximal entangle- 
ment so that this quantum correlation corresponds to 
an ideal electric conductor. When the disorder is intro- 
duced, the entanglement becomes small, the system turns 
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FIG. 2: The distribution of the nearest-neighbor pairwise 
entanglement in the case of weak disorder. The system size 
N = 1600. 
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FIG. 3: The distribution of the nearest-neighbor pairwise en- 
tanglement from weak disorder to strong disorder. The sys- 
tem size N = 1600. 

from extended to localized state, and the conduction of 
system becomes finite and tends to be zero, indicating an 
insulator and no quantum correlation. If we consider the 
thermal conductivity carried by one spin excitation, the 
result is the same. 

In fact, the disorder leads to spatial exponential de- 
crease of eigenstates no matter how weak the disorder 
strength is [20j. When a localization of state occurs, 
what will be the behaviors of pairwise entanglement? 
In the following, we will discuss the entanglement be- 
tween two nearest-neighbor sites, and a distribution of 
the pairwise entanglement is expected. The concurrence 
for the nearest-neighbor sites i and i + 1 is given by 
Ci = 2\ipi?pi + i\, and the numerical results are given in 
Fig. 2. 

In Fig. 2, we notice that the entanglement distributes 
among every pair sites uniformly when the disorder is not 
present. Once the disorder is introduced, the localization 
emerges and the entanglement distribution becomes site- 
dependent. We can see that even the strength of disorder 
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FIG. 4: (a)The maximal pairwise entanglement as a function 
of disorder strength. N = 1600 and 7 = 0.1. (b) The critical 
value of A as a function of relative index j which label the 
distance from one site to the localization center site. 

is very weak, for example A = 0.01, the entanglements be- 
tween most of pairs are suppressed to a very small value 
and only some pair of sites exhibit higher pairwise entan- 
glement around localization center. In other words, the 
entanglement is constrained to some certain pairs of sites. 
The distribution of pairwise entanglement with different 
values of A are plotted in Fig. 3. When A increases, the 
entanglements between more pairs of sites become close 
to zero, and the number of pairs which have enhanced 
pairwise entanglement becomes small. In other words, 
the width of localized entanglement distribution becomes 
small and entanglement is constrained. We can see that 
the localized entanglement is clear in the entanglement 
distribution picture. 

We also notice that the pairwise entanglement between 
localization center site and its near sites are firstly en- 
hanced in the weak disorder regime and then suppressed 
in the strong disorder regime. This property is shown 
in Fig. 4(a), the entanglement Ci i 0+ \ is plotted as the 
function of A, where iq denote the localization center site. 
There exists a maximal value of entanglement at A = 0.8. 
If A < 0.8, the entanglement is increasing. When A is 
beyond 0.8, the entanglement decreases linearly. Such 
behavior of the entanglement can be understood as fol- 
lowing. On one side, the increase of A will lead to entan- 
glement localization, thus enhance Cj D j D -)_i. On the other 
side, from Hamiltonian we know that the increase of 
A will make the term V^cjci predominant, which sup- 
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FIG. 5: The pairwise entanglement between the center site of 
localization and j'-th site, having distance I = \xj — Xi \ from 
center site. Solid line is the fit curve of exponential decay. 
The size of the system N = 1600. 



presses entanglement in the system. The competition be- 
tween the two roles played by the disorder strength leads 
to the maximum value of the concurrence Ci i +i. We 
also study pairwise entanglements between non-nearest- 
neighbor sites and the center sites. They also present 
the maximal behavior during the increasing of disorder 
strength, but the critical value of A will decrease when 
the site is more far from the center site, which is shown 
is Fig. 4(b). A c is plotted as the function of the index j, 
where j is the relative index from the localization center 
and is defined as j = i — io, io is the index of localization 
center. We note that the decrease of critical A is also ex- 
ponential, which coincide with the the exponential decay 
of wavefunction. 

We also consider the pairwise entanglement between 
the center site and it's near sites when the disorder 
strength is fixed. In this case, the concurrence is given by 
Cj = 2 IV'j+io V'io Ij where the label j is the relative index 
from the localization center. As the amplitude of the wave 
function has the maximal absolute value at Xi , by con- 
sidering Eq. (7), Cj becomes 2|V> I 2 1 exp(—\xj +io -x io |/£), 
which measures the entanglement between the localiza- 
tion center and j-th site away from the center. We cal- 
culate the entanglement and get the result with different 
disorder strengths. As expected, the pairwise entangle- 
ment decreases exponentially, agreeing with the simple 
analytical result. This property is shown in Fig. 5. 

From the above results of pairwise entanglement, we 
can find that when disorder is absent, the pairwise en- 
tanglement between every pair of sites is identical and 
independent on site index so that the nonlocality is com- 
pletely shown. For spin system, one spin is entangled 
with the other one far from it as much as the one near 
it. This agrees with the result that in the extended state, 
the system is a perfect conductor for electric current or 
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FIG. 6: Time evolution of the average concurrence of present 
model with different disorder strength. The two extreme ini- 
tial condition is applied (a) Zero entanglement initial state, 
(b) The maximal entanglement state. The system size N = 
1600. 

spin transport. If the strengthen of disorder is not zero, 
the pairwise entanglement becomes localizable and de- 
pendent on site index. The entanglement between near- 
est neighbor spins (sites) is different, and spins (sites) will 
prefer to entangle with near spins (sites) than the one far 
away. The nonlocality is destroyed, and entanglement is 
decrease globally, although some local entanglements be- 
tween pair spins do not. In this case, one spin excitation 
would be pinned at some sites and the conduction van- 
ishes. 

We have studied the properties of ground-state entan- 
glement in the one-dimensional Anderson disorder model. 
Next, we will investigate the dynamics of entanglement. 
This question can be solved by calculating the time- 
dependent Schrodinger equation 

= - i/H-i + V^i (10) 

which can be integrated numerically. 

The dynamical behavior of average concurrence of pre- 
senting model are shown in Fig. 6 with different disor- 
der strengths. Fig. 6 (a) displays the dynamical evolu- 
tion of entanglement with the initial state being \ip)o = 
|1) = cj|0), where I, indicating one site of system, is set 
I = N/2 and N is the size of system. In this case, the 
initial state has no entanglement, what will this state 
get entangled when time goes on? If disorder is ab- 
sent, A = 0, the average concurrence increases linearly 
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with time. When time is long enough, the entanglement 
reaches a maximal value, after that it will decrease to a 
smaller value. When the disorder strength is not large, 
for example, A = 0.05, the entanglement also increase lin- 
early at first stage and reach a maximal value. But the 
decrease of entanglement after maximal value becomes 
small. We notice that when A = 0.1, the entanglement 
does not decrease, and almost keeps the maximal value. 
When disorder strength is large enough, like A = 0.5, the 
initial entanglement increases slowly. After reaching a 
small value, it does not increase again. If A is larger, for 
example (A = 1.0) , the entanglement exhibits almost no 
increase. This is a strong indicator of localization. 

If one start the Equ. (jlOJI from the maximally pairwise 
entangled initial state, namely the W state |2|,|2l|, the dy- 
namical behaviors are displayed in Fig 6. (b). The max- 
imal entanglement is invariant while disorder is absent. 
But as soon as the disorder is introduced, the entangle- 
ment decreases as time goes on, the rate of the decrease is 
determined by the disorder strength. The larger disorder 
strength will make entanglement decrease quicker than 
smaller one. But it is interesting that no matter how 
large the disorder strength is, the entanglements will all 
decrease to a same residue value. Comparing the dynam- 
ical evolution of the entanglements from two different ini- 
tial states, we find the asymptotic behaviors are different. 
The residue values of initial maximal entangled state are 
larger than those of initial non-entangled state for any 
value of disorder strength. 

The evolution of entanglement in this system is deter- 
mined by the diffusion of wavefunction and the localiza- 
tion. Without disorder, the wavefunction of electron will 
spread along the spatial direction. The external disorder 
potential can localize the motion of electron (or one exci- 
tation in XY spin system) on finite region of space. The 
existence of these two process leads to the dependence 
of evolution of entanglement on the initial condition and 
disorder strength. For an initial unentangled state which 
has a strong tendency to spread, a weak disorder can not 
prevent occupation possibility of particle sharing among 
more sites, so the entanglement will increase a lot, but 
a strong disorder will make a localization of the state, 
which leads to a small increase of entanglement. If the 
initial state is maximally entangled state, a completely 
extended state, so only the localization will determine 
the evolution entanglement. The entanglement must de- 
crease as time goes on and the disorder strength deter- 
mines the rate of decrease. We also do same calculations 
when we use other initial states different from the above 
two states. The result is also the same, namely, the disor- 
der will destroy entanglement. By studying these prop- 
erties, we can find that other states, except for ground 
state, have similar entanglement characters. So, it is true 
that electric conduction with single electron approxima- 
tion and thermal conduction of one spin excitation are 
connected with quantum correlation measured by pair- 



wise entanglement. 

In summary, we have studied the ground state and 
dynamical pairwise entanglement of one-dimensional An- 
derson disorder model. By simple analytic and numerical 
results, we can find the entanglement is the indicator of 
localization caused by disorder. The disorder can destroy 
such quantum correlation. The relations between the en- 
tanglement and conduction are also discussed. On the 
other hand, we can localize one qubit on certain site by 
disorder, then we can do quantum operation on it. 

It is interesting to consider entanglement in other dis- 
order models. These studies will strength our under- 
standing of entanglement, localization, and their rela- 
tions. For instance, it is an interesting question to study 
effects of disorder on entanglement, the relation between 
localization and entanglement for other subspace with 
the number of electrons or spin excitations being large 
than one, which are under consideration. In that model 
, the effect of disorder on multi-party entanglement can 
also be investigated. 
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